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Abstract 

We study the 4(i Heisenberg spin glass model with Gaussian nearest- 
neighbor interactions. We use finite size scaling to analyze the data. We 
find a behavior consistent with a finite temperature spin glass transition. Our 
estimates for the critical exponents agree with the results from e-expansion. 
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The lower critical dimension di of the short range models remains one of the most 
controversial questions in the spin glass theory. In the last ten years a consis- 

tent number of works has been devoted to studying the Ising model, which seems to 
be 1^ very close to di in d = 3. In the case of the short range isotropic Heisenberg 
spin glass, the conclusions of various computer simulations in d = 3 0] agree that 
the system is below di. Using domain wall renormalization group techniques, it was 
argued [|, || that di = 4 for this model. To our knowledge, there are no previous 
numerical simulations on the 4d Heisenberg spin glass apart from an early work by 
Stauffer and Binder [§], which studied the time dependence of the Edward Ander- 
son order parameter qEAif) for vector spin glasses from d = 2 to d = Q. They 
observed a change in the behavior of the Heisenberg model in d = A but concluded 
that probably there was no finite T phase transition. 

We have studied the Ad isotropic Heisenberg spin glass model with Gaussian 
nearest-neighbor interactions. We have simulated small lattices (with linear size 
from L = 3 to L = 5), using finite size scaling to analyze the data. Our main 
conclusion is that d = 4 seems to be well above the di of the model. 

After defining the model and the quantities we have measured, we present a 
theoretical estimate for critical exponents obtained from e-expansion results [|l2 



and we discuss the computer simulations in some details. Finally, we present the 
numerical results and the conclusions. 

The Hamiltonian of the model is given by 

1 ^ 

where the spins {(fj} belong to the unit three-dimensional sphere. The sum runs 
over the nearest-neighbor pairs in a simple hypercubic lattice with periodic boundary 
conditions and size V = L'^. The interactions Jij are quenched independent random 
variables with a symmetric Gaussian distribution 

P{J,,) = (2) 

The spin glass correlation function can be defined as 

GsGir)^^ E (3) 

\ri-rj\=r 

where < (■) > means thermodynamic average and (■) means average over samples. 



1 



In the thermodynamic hmit, in the paramagnetic phase, for distances r >> ^sg 
we expect Gsci^) oc e~^^^'^'^, where we have introduced the spin glass correlation 
length ^sG- The divergence of ^sg when the critical temperature T^. is approached, 

^SG^iT-T,)-^ (4) 

is characterized by the critical exponent u. The power law decay of the correlation 
function for large r at the critical point, 

GsGir) oc ^ (5) 

is described by the anomalous dimension r]. 

In order to understand the behavior of the model, it is interesting to introduce 
two independent replicas of the system, with the same disorder configuration. Since 
we are dealing with a vector spin glass, the overlap is a second-rank tensor in spin 
space. We can define 

Q'-'^^E^^ri: (6) 
^ 1=1 

where {a} and {r} are the spins of the two replicas and /i, z/ = 1, 3 refer to spin 
components. In our case the system is completely isotropic, therefore we expect 
to deal with one spin glass order parameter q. It is convenient to consider the 
rotational invariant quantity 



1 3 



. TT E il"')' (7) 



Q 

The order parameter probability distribution P{q) is correspondingly given by 



Piq)^<6iq-Q)> (8) 
One of the quantities we have calculated is the spin glass susceptibility, defined 



as 




XsGiL,T) = -< [l^cj.n] > (9) 



It can be easily verified that this definition is equivalent to xsg = < >, the 
factor 3 having been inserted to obtain limT^oo Xsg{.L, T) = 1, like in the Ising case. 
In the thermodynamic limit we expect the spin glass susceptibility to diverge when 
Tc is approached from the paramagnetic phase as 

X5g(T)oc(T-T,)-^ (10) 
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where, by hyperscaling, 7 = (2 — r^) i/. 

If Tc = and the system is below di, we still expect C,sg(T) and xsg(T) to 
behave near the critical point according to (^) and ( ]T(]| ) respectively, diverging with 
power laws. We expect exponential divergences for d ^ di, because i/, 7 00 in 
this limit. In the 3d Ising spin glass the difficulty in distinguishing between the 
system being at or above di |Q is correlated to the large value of u obtained from 
finite size scaling analysis. We will see that this seems not to be the case in our 
model. 

The Binder parameter has proved very successful to establish the presence or the 
absence of a finite T phase transition. Its definition for the Ising spin glass model 



njl can be easily extended to the Heisenberg case. In the high temperature region, 
where we can neglect interactions, the q'^'^ are approximately independent variables 
with the same symmetric Gaussian distribution of width ~ V^^^"^. The function 



g(L,T) = -(n-9 ^^'^ -) (11) 

is a dimensionless parameter defined so that g < I. From an explicit calculation, we 
have obtained g ~ 1/V for T — * 00. In the thermodynamic limit we expect therefore 
g(T) = above Tc and g(T) = 1 for T = 0. If there is a T = singularity we expect 
the curves of g{L, T) against T for different L to come together as T — > 0, while for 
a finite T transition we expect them to intersect at Tc, which allows to locate the 
critical point quite precisely. 

The mean field spin glass critical temperature T^^^ for m-component spins be- 
longing to the unit sphere, with coordination number z and Jfj = 1, is approximately, 
for large z, T^^ ~ ^jm. In our case m = 3 and z = 8, so 

T^^ ~ 0.94 (12) 

This value is a reasonable upper limit to the temperature range in which we can 
expect a phase transition. 

If there is a nonzero Tc, the behavior of the system at the transition point is 
characterized by two independent critical exponents. We have obtained a theoret- 
ical estimate for the spin glass correlation length exponent v and the anomalous 
dimension 77 from e-expansion results. This is possible in the Heisenberg case be- 
cause the coefficients, calculated to the third order by Green |jl^ , show the expected 
oscillatory behavior: 



ri = -0.2 £ + 7. 7333- 10^2 ^2 _ 78127. 10-2^3 + 0(e^) 
+ = -1.2 £ + 1.164 1.4735 e3^0(£^) 
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V 


ri 


0.7096 


0.8657 


-0.2256 


-0.4945 


0.8224 




-0.2976 





Table 1: v and 77 from the corresponding V\N N being the row and D the 
column index. 



where e = du — d = Q — d, du = 6 being the upper critical dimension of short range 
spin glass models. 

We have used the simple resummation method of the Pade approximants, in 
which the series is replaced by the ratio V[N, D] of a polynomial of degree to one 
of degree D. The values obtained for e = 2 are presented in [Tab. 1]. By comparing 
results from V[l, 2] and V[2, 1] we can estimate, for the 4c? Heisenberg spin glass, 

z/ ~ 0.8 r/ ~ -0.4 (14) 

with a larger uncertainty on the value of 77. 

We have simulated hypercubic lattices in 4d with periodic boundary conditions 
and linear sizes L = 3, 4 and 5. The number of samples is 400, 200 and 100, 
respectively. Simulations have been performed in the region T < 1, down to T = 0.4 
for L = 3, to T = 0.45 for L = 4 and to T = 0.5 for L = 5. We will see that T = 0.5 
seems to be very close to the of our model. 

In order to thermalize our samples, we have used Simulated Tempering |ll2|-||14||. 



already proved very efficient for Ising spin glasses 0, The system, in our case 
consisting of two independent replicas, is allowed to change temperature between a 
fixed set of {/5„}, where we can take for simphcity Pn+i > Pn, P = 1/7" becoming a 
dynamical variable. The stationary probability distribution for the configuration C 
at Pn is given by 

Poo(C, n) oc e-^tot[C, n] ^^^^[^^ ^ ^^^^^j _ (^5) 

where we have defined the extended Hamiltonian T-Ctot[C,n], {gn} being a set of 
arbitrary numbers chosen a priori. After reaching the equilibrium, the system 
moves between the remaining at the equilibrium. In order to obtain the same 
stationary probability for all the different temperatures we have to take Qn = PnFn, 
Fn being the total free energy at This means that, for n' = n ± 1, at the first 
order in A/3 = {Pn' - Pn)- 

AHtot = Apn - {gn' - gn) ^Ap(n- i(< Hn' > + < Hn >)) (16) 
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where Ti. is the instantaneous value and < Tin > is the statistical expectation value 
at j3n of the total energy. 

Our Simulated Tempering steps are just Monte Carlo steps at the end of which 
the system is allowed to change temperature, the new /5„ suggested being Pn±i with 
equal probability. Obviously during the MC step the two rephcas evolve indepen- 
dently. Since we are dealing with a model with continuous degrees of freedom, there 
is an arbitrary parameter in the Metropolis algorithm, corresponding to the maxi- 
mum rotation angle 9max permitted to single spins in one step. We have chosen 9max 
in order to obtain the acceptance as close to 1/2 as possible. 

We have been careful in fixing the set of temperatures for the different sizes. 
Two contiguous values of have to be as different as possible to help in decorre- 
lating without making too small the corresponding transition probability. We have 
used as a basic criterion the condition that there was a non-negligible overlap in 
the values of the energy computed at contiguous /?„ for each sample. In our case 
this was verified by choosing equidistant temperatures, with AT = 0.1 for L = 3 
and AT = 0.05 for L = 4 and 5. In order to perform simulations down to lower 
temperatures, particularly in the L = 5 case, it would be necessary to decrease AT 
with a considerable greater amount of computer time. 

We have used a slow cooling procedure to take the system near the equilibrium 
at the lowest temperature. Statistics were collected over the last part of the about 
3000 L MC steps at each temperature, to evaluate approximately the corre- 
sponding <nn>. The next about 70000 for L = 3, 200000 for L = 4 and 400000 
for L = 5 Simulated Tempering steps were used to thermalize the system and to 



improve iteratively |15] the estimates for the {gn' — gn}- Finally, all the quantities 
we were interested in were calculated in the last part of the Simulated Tempering 
cycle, of about 140000 steps for L = 3, half a million steps for L = 4 and more than 
a million steps for L = 5. 

In Simulated Tempering, we can estimate the statistical expectation value 
< On > of an observable O at /5„ as 

< >= C'(t)5/3(J)/3„ fn = Ji=Yl ^m^n (17) 

where we have defined the frequency /„ observed at Pn, being the total number of 
steps. With the given choice for the {gn}, we expect /„ ~ l/rir, tit being the total 
number of temperatures considered (in our case, = 7 for L = 3 and = 8 for 
L = 4 and 5). The closeness of the /„ observed at different temperatures represents 
one of the main verifications that Simulated Tempering works well. In the last part 
of the cycle we have obtained max{/„} < 4min{/„} for each sample, the /„ being 
compatible with the expected values for the different ut- We have also checked 
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how many rie times the system was moving from one extreme to the other of the 
temperature range, obtaining for each sample rie > 100 in the L = 3 case and 
fig > 200 for L = 4 and 5, reasonably large values for the system really exploring 
the entire phase space. 

In order to check thermalization, besides verifying that the < q'^" > were com- 
patible with zero for each sample, we have divided the last part of the Simulated 
Tempering cycle in five equal intervals, checking that the various computed quan- 
tities show no evident drifts. This was verified also for the spin glass susceptibility 
(H) and the fourth moment of the P{q) (||) in the L = 5 case. 

The simulations have taken in all about 6 months of DecSOOO-workstation. 

Our numerical results for the spin glass susceptibility @ and the Binder pa- 
rameter (|TT|) are presented in [Fig. 1] and [Fig. 2] respectively. The behavior of the 
curves of g{L,T) as a function of T for the different lattice sizes strongly suggests 
the presence of a finite T spin glass transition, the noncoincidence of the intersection 
points being presumably due to systematic corrections to finite size scaling. In order 
to confirm that curves really splay out below the critical point, as expected in the 
presence of long range spin glass order, it would be obviously preferable to obtain 
data on the g{L,T) down to lower temperatures and for larger lattice sizes. 

If scaling is satisfied, near Tc we expect 

Xsg{L,T) = L'-^xsg{{T-T,)L'/'') (18) 
gi^T) = ~g({T-T,)L'lA (19) 



where xsg and g are scaling functions while v and t] are respectively the spin glass 
correlation length exponent (^ and the anomalous dimension (|^) previously defined. 

From the xsc scaling law (p!8D, using a standard three-parameter fitting routine, 
we have obtained 

Tc = 0.50 ±0.06 z/ = 0.61 ±0.08 2 - r/ = 1.8 ± 0.5 (20) 
while requiring that the Binder parameter data scale according to equation (|19|), 



with a more simple two-parameter fit, we have found 

Tc = 0.52 ±0.02 i/ = 0.89 ±0.06 (21) 

In [Fig. 3] and [Fig. 4] we present the corrispondent scaling plots. It must be 
emphasized that statistical errors quoted here are just a delimitation of the range of 
values beyond which our data do not scale well. Systematic errors due to corrections 
to finite size scaling cannot easily be evaluated but could be quite important, because 
of the small lattice sizes considered. 
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Our results agree well with a nonzero Tc for the Ad short range Heisenberg spin 
glass (|l]), giving in this case of Gaussian nearest-neighbor interactions @ the value 
~ 0.5, that is however well below the T*^^ ~ 0.94 (^) of the model. 

Systematic corrections to finite size scaling may explain the discrepancy between 
the estimates for u obtained from the xsg ( pO|) and the Binder parameter (|2l|) 



respectively. Our most significant result is nevertheless that u seems to be not 
large, being both values smaller than 1. As we have already pointed out, a large 
value of u would suggest the system being at ~ c?;, since we espect i/, 7 — > oo for 
d ^ di. In this case, therefore, the result z/ < 1 is consistent with the short range 
Heisenberg spin glass model being well above di m. d = 4. 

Comparing the two estimates for i/, we can obtain the approximate value 
V ~ 0.75, that agrees well with the theoretical estimate (|I^) v ~ 0.8. From the 
hyperscaling law a = 2 — z/ we can also estimate the corresponding value of a, 
describing the critical behavior of the specific heat. We find a ~ — 1, accidentally 
not far from the mean field value. 

Finally, our statistics are inadequate to obtain a significant estimate for rj. The 
value found from the xsg data (pOf), 2 — 17 = 1.8 ± 0.5, is however compatible with 
the theoretical estimate ( p^lf ) 7] ~ —0.4, even if not in good agreement. 

We have simulated small lattices, with a rather small amount of computer time. 
More accurate simulations are necessary in order to confirm the presence of a finite 
T phase transition in the Ad short range Heisenberg spin glass model, improving our 
estimates for critical exponent. Simulated Tempering seems to be highly suitable for 
this purpose. It might be also interesting to devote some attention to the behavior 
of the to our knowledge not yet extensively studied in the case of short range 

vector models. 

I am deeply indebted to G. Parisi for his continuous support. I would also like 
to thank E. Marinari for many helpful discussions, A. Baldassarri, F. Ritort and 
J. J. Ruiz for their useful suggestions. 
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1 1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 

Figure 1: The spin glass susceptibility xsg as a function of (5 for the different lattice 
sizes. Lines are only to join neighboring points. 




Figure 2: The Binder parameter 51 as a function of T for the different lattice sizes. 
Lines are only to join neighboring points. 
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0.5 1 1.5 2 2.5 
(T - 0.50)LVO-62 

Figure 3: The scaled spin glass susceptibility xsg/ L'^~^ as a function of the scaled 
reduced temperature (T — Tc)L^^'' . The line is just guide to the eye. 




-0.5 0.5 1 1.5 

(T - 0.52)LVO-89 



Figure 4: The Binder parameter as a function of the scaled reduced temperature 
(T - Tc)L^'''. The line is just guide to the eye. 
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